Anisotropy of in-plane magnetization due to nodal gap structure in the vortex state 
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We examine the interplay between anisotropy of the in-plane magnetization and the nodal gap 
structure on the basis of the approximate analytic solution in the quasiclassical formalism. We 
show that a four-fold oscillation appears in the magnetization, and its amplitude changes sign at 
an intermediate field. The high-field oscillation originates from the anisotropy of the upper critical 
field, while the low-field behavior can be understood by the thermally activated quasiparticles near 
nodes depending on the applied field angles. The temperature dependence of the magnetization also 
shows a similar sign change. The anisotropy of the magnetization ofi'ers a possible measurement to 
identify the gap structure directly for a wide class of type II superconductors. 
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In the last two decades, a family of unconventional su- 
perconductors has extensively grown up, and it has cre- 
ated more demand for a direct observation of the gap 
structure. Although thermodynamic quantities at low 
temperature follow power-law behaviors due to the pres- 
ence of gap nodesUQ, this in principle cannot determine 
an absolute direction of nodes, which has a significant im- 
portance to selectively identify the pairing mechanism. 
Recently, a powerful method has been developed to mea- 
sure the gap nodes directly. Namely, field angle depen- 
dences in oriented magnetic fields have shown oscillatory 
behaviors in the thermal conduct ivitvH, 0, IE la| 
and the specific heat d El m El El 

Il4l | , reflecting low- 
energy excitations inherent from the gap structure. 

An analysis of those two experiments, however, is 
somewhat ambiguous, and leads opposite conclusions at 
worstHEl- This is because a dominant source of the 
oscillation differs depending on what fields and tempera- 
tures we consider. For instance, at low temperature the 
activated quasiparticle (QP) due to the supercurrent en- 
ergy shift roughly determines the oscillatory behaviors in 
the specific heat |15tll6t.ll7i . On the other hand, a trans- 
port lifetime in the thermal conductivity plays an impor- 
tant role in intermediate temperature range, which could 
reverse the oscillation caused by the density of states 
(D0S)E1E1|23- Therefore, more such angle-resolved 
bulk measurement is highly required to draw a definite 
conclusion with a consistency among different probes. 

In this Letter, we examine the interplay between 
anisotropy of the in-plane magnetization and the nodal 
gap structure in the whole region of the H-T diagram, 
and discuss whether the magnetization is a possible probe 
to identify the gap structure or not. Experimentally, the 
basal plane anisotropy of magnetization, M{(f)h), has al- 
ready been observed for example in non-magnetic boro- 
carbides, RNisBzC (R=Lu, Y)0, i^. A remarkable 
fact is that the amplitude of the oscillation in M changes 
sign from high to low fields [2ll|. Temperature dependence 
also exhibits a similar sign change |22|. It is worth noting 
that the low-_ff and low-T signs are opposite to what is 



expected from the known angular dependence of the up- 
per critical field, Hc2- Indeed, such a sign change has a 
close resemblance to that observed in the specific heat of 
Sr2Ru04|l Q, where the sign change can be ascribed 
to the modulation of the gap amplitudeflTj. Hence it 
should be squarely addressed. 

In the tetragonal symmetry the standard London the- 
ory does not introduce any in-plane anisotropy to the 
mass tensor. Kogan et al. argued that the nonlocal cor- 
rection to the London theory gives rise to anisotropy in 
the mass tensor, and qualitatively reproduced the overall 
observed behaviors |22j. Nevertheless, it is natural to ad- 
dress that the nodal gap structure alone could create the 
anisotropic behavior in Al. Theoretically, the validity of 
the nonlocal London theory is also questionable deep in 
the vortex state, where an overlap of the cores plays a 
dominant role. In the high-field end the anisotropy of 
Hc2 predominates over any other contributions. 

Our approach is based on the approximate analytic so- 
lution of the quasiclassical equations near Jf^^lisl li^ . 
In this approximation, we replace the spatial depen- 
dence of the magnetic induction, B{R) and the diago- 
nal quasiclassical propagator, gk{R) = — * / dS^G/n by 
their average, B and gk, and the Abrikosov solution, 
A(fc, R) is used for the vortex lattice structure. We 
consider the two-dimensional cylindrical isotropic Fermi 
surface for simplicity, and focus on the (ia,2_j,2-wave, i.e., 
A(fc) — Aip{k) with ip{k) = •\/2cos(2</)), where cj) is the 
azimuthal angle of k measured from the a;-axis. Then, 
we have the analytic solution of the propagator as 
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where u„ — 2uJn/v±{k)y^2\e\B with the fermionic Mat- 
subara frequency and W{z) ~ e^^^erfc(— iz) is the Fad- 
deeva function. Here v± is the component of the Fermi 
velocity perpendicular to the field, and v'j_ (x 1-1-2 sm^{(f>— 
(/i/i — 7r/4) for the in-plane field measured from the nodal 
direction. The prefactor of the velocity characterizing 
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the quasi two-dimensionality can be absorbed into the 
definition of Hc2- Note that only the spatial average of 
the gap, A(fc) = A(fc, R) appears in the expression. The 
solution recovers the correct BCS result in the limit of 
B 0. 

Although the approximation seems to be valid near 
Hc2: the present author has shown that the analytic so- 
lution captures an essential physics in low-iJ and low-T 
regions as well if we minimize the free energy with re- 
spect to B and A^,!!^,!!^- In the clean limit, the free 
energy measured from the normal state is given bv[2.5) 



lAl^ln (4e'^nc) 
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with 
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V^|A(fc)|2iy(iM„), (3) 



where A^'o is the DOS in the normal state, 7 ~ 0.5772 is 
the Euler's constant and (• • • ) denotes the average over 
the Fermi surface. The cut-off frequency = lUc/'^t^T = 
60 is used in the present paper. Note that the dimen- 
sionless quantity, H^2/-^o'^c ^ proportional to (A/^)^, 
so that we define it as our Ginzburg-Landau (GL) pa- 
rameter . 

Figure ^ shows the magnetization curves for the field 
parallel to the nodal direction, i.e., (ph = 0. The GL 
parameter is k/kc — 5, where Kc — 5.363 is the criti- 
cal value between type I and II superconductors. The 
calculated curves well describe the correct iJ-linear de- 
pendence near Hc2 and the logarithmic increase at low 
fields. A smooth peak rather than the cusp-like behavior 
near Hd is obtained due to the artifact of the approxima- 
tion. The inset shows the H-T phase diagram, where the 
Hci is determined by the field showing the maximum in 
the magnetization curve. The Meissner state is realized 
below Hci, namely, the free energy has a minimum at 
B = 0. The theory well describes the low-field behavior 
at low temperatures, which is the most severe region for 
the present approximation. 

Next we discuss the field-angle dependence of the mag- 
netization. Figure 121 shows — M as a function of (ph for 
several H at T/Tc = 0.1. It is remarkable that the 
nodal gap alone gives rise to the four-fold oscillation of 
—M{(j)h), which is reversed as H decreases. Note that 
the data at different H are vertically shifted to fit in a 
single figure. The inset shows the anisotropy of Hc2i4'h), 
which is consistent with the known results. In the present 
approximation at T = and H < Hc2 , we have 




FIG. 1: The magnetization curves for several temperatures 
with k/kc = 5, where Kc is the critical value of the GL pa- 
rameter K between type I and II superconductors. The inset 
shows H-T diagram, where Hci (T) is determined by the field 
showing the maximum in the magnetization curve. 
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FIG. 2: The magnetization as a function of the field angle 
measured from the nodal direction. The data at different H 
are vertically shifted. The inset shows oscillation in Hc2{4>h)- 
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where A{<j)h) = {\>p{k)\'^ / v\{k)) > 0. From these expres- 
sions, we expect the similar oscillatory behaviors between 
—M{(f)h) and Hc2{4>h) at high fields. 

To quantify the amplitude of oscillations, we define 
the amplitude as -5M = -[M{(j)h = 7r/4) - M{(j)h = 0)]. 
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FIG. 3: The amplitude -5M = -[M((;/>h = 7r/4)-Af(<;/>h = 0)] 
of the four-fold oscillations of the in-plane magnetization as 
a function of (a) H at T/Tc = 0.1 and (b) T at HjHc^ = 0.4. 
The amplitudes of the spatially averaged gap magnitude, (5A, 
and the number of the activated quasiparticles near nodes (see 
text in detail), (Jrinode, are also shown. 



The field dependence oi—SM for T/Tc — 0.1 is shown in 
Fig.Ofa). The amplitude —SM decreases almost linearly 
as H decreases and changes its sign a.t H = H* . It has 
already shown that the Hc2 anisotropy brings about the 
positive amplitude in high fields. The sign change in de- 
creasing field indicates that different physics comes into 
play to predominate over the effect of the Hc2 anisotropy. 

To elucidate this point, the H dependence of the os- 
cillatory amplitude of the gap magnitude, (5A and the 
population of the thermally activated QPs near nodes, 
5n^iodc are shown in Fig. Ofa). The amplitudes are de- 
fined in the same manner as —SM. We have defined nnodc 
normalized by the normal-state counterpart as 
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j;;^dLofiLu/T)N{uj,k,±s) 
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where N{lu, k) ^ Ke g is the DOS in the superconduct- 
ing state, f{x) is the Fermi-Dirac distribution function, 
and k, = (7r/2)i + 7r/4 (i = 0,1,2,3) with 6 = 0.027r. 
Note that the field-angle dependence is implicitly in- 
cluded in N{uj, k). Since the gap amplitude has a strong 
H dependence near Hc2, i.e., ^{4>h) oc \Jl — H/Hc2{(t>h), 
the anisotropy in Hc2{4'h) yields the larger anisotropy in 




FIG. 4: The temperature dependence of the crossover field 
of reversible oscillation in the magnetization for several GL 
parameters. 



A(0/i) and hence the larger —M{(j)h) at closer Hc2 as 
shown in Fig.l^ta). Note that the relation between —M 
and A is given by -M(0ft) = NqA^ {(/),,) /2Hc2i(l)h) at 
T = near Hc2- The anisotropy in A becomes less pro- 
nounced as H decreases. 

In contrast to the H dependence of 6 A, the positive 
^n-nodc increases as H decreases. In other words, more 
quasiparticles are thermally activated near nodes when 
the field is applied in the anti-nodal direction as com- 
pared with the case of || fcnodc. This is understood by 
the so-called Doppler shift argument as follows, which is 
valid for a low density of vortices at low temperatures 2C|. 
The local quasiparticle with the momentum k has the 
energy spectrum + Vg ■ k in the supercurrent fiow- 
ing around vortices with the velocity Vg. This energy 
shift gives rise to a finite DOS near the gap nodes if 
\k ■ Vs\ > |A(fc)|. Since -L H, the number of the rele- 
vant gap nodes becomes largest in the case of i? _L fcnodc- 
Thus, we have Sunodc > 0. Then, the reduction of the su- 
perfluid density leads to increase of the penetration depth 
as well as the magnetic induction. As a result, we have 



smaller —M for H _L k„ 



idc 



than that for if 11 fen 



The temperature dependence of — M, ^A and 5n,iodc 
shows a similar behavior as shown in Fig. |3Jb), hence, 
a similar argument holds for T dependence as well. We 
then conclude that in the presence of the gap nodes, a 
sign change should occur in the amplitude of — M as a 
function of H and T. The sign of the oscillatory ampli- 
tude has an important information about the direction 
of gap nodes. 

Finally, we discuss the dependence of the GL param- 
eter, K. Figure 0] shows the temperature dependence of 
the crossover field, H*, for several k's. The crossover 
field is located in the intermediate field range as a con- 
sequence of the competition described in the above. It 
is almost independent of n at least for k/kc > 5. This 
implies that the identification of the gap structure us- 
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ing the field-angle dependence of the magnetization can 
be applied for a wide class of the intermediate and the 
strong type II superconductors. 

In summary, we have shown the interplay between the 
in-plane anisotropy of the magnetization and the nodal 
gap structure with four line nodes. The field-angle de- 
pendence of the magnetization, M, exhibits a clear four- 
fold oscillation, as a direct consequence of the nodal gap 
structure. The high-field oscillation in —M shows a min- 
imum when the field is applied along the nodal direc- 
tion, while the low-field one has a maximum in the same 
field direction. The sign change in the oscillatory am- 
plitude can be understood by the competition between 
the effect of the Hci anisotropy and the low-field quasi- 
particle excitations near nodes at low temperatures. As 
a result, the crossover field appears in the intermediate 
field range, and is almost independent of the Ginzburg- 
Landau parameters. The temperature dependence shows 
a similar sign change. All these features suggest that the 
anisotropy of the magnetization can be a good candidate 
to identify the gap structure directly for a wide class of 
type II superconductors. 
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